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Abstract. A Filbert matrix is a matrix whose entry is l/Fi+j—i, where 
F n is the n th Fibonacci number. The inverse of the n X n Filbert matrix 
resembles the inverse of the n X n Hilbert matrix, and we prove that it shares 
the property of having integer entries. We prove that the matrix formed by 
replacing the Fibonacci numbers with the Fibonacci polynomials has entries 
which are integer polynomials. We also prove that certain Hankel matrices 
of reciprocals of binomial coefficients have integer entries, and we conjecture 
that the corresponding matrices based on Fibonomial coefficients have integer 
entries. Our method is to give explicit formulae for the inverses. 



1. Introduction 

The n x n Hilbert matrix is the n x n matrix whose (i,j)-entry is i+ j_ 1 . In 
H|, Man-Duen Choi explores many fascinating properties of the Hilbert matrix, 
including the fact the the (i, j)-entry of its inverse is 

Choi asks what sort of coincidence it is if the inverse of a matrix of reciprocals of 
integers has integer entries. In this paper we show that the inverses of the Hankel 
matrices based on the reciprocals of the Fibonacci numbers, the reciprocals of the 
binomial coefficients , and the reciprocals of the binomial coeffiencts ( I+ 3 +2 ) 
all have integer entries. We also find formulas for the entries of the inverses of these 
matrices and related matrices. 

Definition 1.1. Let {a^} be an integer sequence with a& ^ for A; > 1. A recip- 
rocal Hankel matrix based on {dfc} is a matrix whose (i,j)-entry is l/aj+j_i,. We 
denote the n x n reciprocal Hankel matrix based on {a&} by R n (ak)- 

The formula for the entries of the inverse of R n {Fk) bears a striking resemblence 
to the formula for the entries of the inverse of the Hilbert matrix. Therefore, we 
call a reciprocal Hankel matrix based on the Fibonacci numbers a Filbert matrix. 

2. Filbert matrices 

We need the Fibonomial coefficents to describe the inverse of the Filbert matrix. 
See [til for more information on the Fibonomial coefficients. 



Definition 2.1. The Fibonomial 
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where n and k are nonnegative integers. 

Theorem 2.1. Let e(n,i,j) = n(i +j + 1) + (|) + Q + 1, and Zet &e the 

nxn matrix whose (z, j)-entry is 

';":ir))((": j :; 1 ))(Ci-T 2 )) 2 



^■H = (-i) e(n ^' ) ^ +J -il 



T/ien the nxn matrix W{n) is the inverse of the Filbert matrix R n (Fk), and W(n) 
is an integer matrix. 



This theorem is a special case of Theorem 2.2, which we prove below. The 
formula for the entries of the inverse closely corresponds to the formula for the 
entries of the inverse of the nxn Hilbert matrix. It results from Eq. (|l]) by changing 
all binomial coefhecients to Fibonomial coefficients and changing the exponent of 
— 1. The pattern of the signs of entries the inverse of R n {Fk) is that they are 
constant on 2 x 2 blocks, and alternate between blocks. 

The Fibonacci polynomials f n {x) are defined by fo(x) — 0, fi(x) = 1, and 
fn{x) = xf n -i(x) + f n ~2(x) for n > 2. We also use /„ to denote the Fibonacci 
polynomial f n (x), especially when we want to reduce the clutter in some equations. 
The JE-Fibonomial coefficients are the obvious generalization of the Fibonomial 
coefficients. 

Definition 2.2. The x-Fibonomial coefficients are 

fn-i+l{x) 




n 



i=i 

where n and k are nonnegative integers. 



fi(x) 



To form the (i, j)-entry of the inverse of R n {fk{x)), replace each Fibonacci num- 
ber and Fibonomial coefficient in Wij (n) with the corresponding Fibonacci poly- 
nomial and ir-Fibonomial coefficient. 

Theorem 2.2. Let V(n) be the nxn matrix whose (i, j)-entry is 




Then the nxn matrix V(n) is the inverse of the Filbert matrix R n (fk{x)), and the 
entries of V(n) are integer polynomials. 

The recurrence ((£))= /n^ff 1 ^)) + fn-k+i(x) ((£)) shows that the 
Fibonomial coeffcients are integer polynomials, which implies that the entries of 
V(n) are integer polynomials. 

3. Technology 



The proof of Theorem 2.2 and proofs of succeeding theorems amount to proving 
various identities involving sums of products of Fibonomial coefficients and binomial 
coefficients. We supply computer proofs of these identities. In some cases, the 
computer cannot do the entire proof directly, and human intervention is required 
to separate the proof into smaller pieces that can be done by computer. 
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The programs and packages used to produce the proofs for this paper include 
Maple V Release 5, the Maple package EKHAD written by Doron Zeilberger, and 
the Mathematica package MultiSum written by Kurt Wegschaider. EKHAD is de- 
scribed in pL and it is available at www. math, temple, edu/^zeilberg. MultiSum 
is available through the web site www.risc.uni-linz.ac.at/software/, and it is 
described in Q. The particular functions that we use from these packages are zeil 
from EKHAD and FindRecurrence from MultiSum. 

Both of these functions find a telescoped recurrence for a summand F(n,k), 
where k is the summation variable. The function zeil uses Zeilberger's algorithm 
to find a rational function R(n,k) and a recurrence operator P(n,N), where N is 
the shift operator in n, such that 

P(n, N)(F(n, k)) = R(n, k + l)F(n, k + 1) - R(n, k)F(n, k). (2) 

Let f(n) be the unrestricted sum F(n,k). In many situations, Eq. (^|) implies 
that P(n, N)f(n) = 0, making it easy to verify that f(n) is constant. 

The function FindRecurrence gives similar results with summands of the form 
F(n, k) where n and k are vectors. 

Maple V Release 5 also includes an implementation of Zeilberger's algorithm as 
the function sumrecursion of the package sumtools. However, sumrecursion only 
gives the recurrence operator P(n, N), and not the rational function R(n, k), which 
will be essential when we prove identities invol ving a restricted sum. 

The sums involved in the proof of Theorem |2.2| are of products of Fibonomials, 
not binomials, so these procedures do not apply. However, we obtained recurrences 
for sums of products of Fibonomials by modifying recurrences found by these pro- 
cedures for the corresponding sums of products of binomials. 

4. Proof of Theorem |2~2] 
The (z,m)-entry of the product V(n)R n (fk{x)) is 

n 

p(n,i,m) = ^2P(n,i,m,j), 
i=i 

where 
P(n,i,m,j) = 

The summand satisfies the following recurrence relation that is related to a re- 
currence produced by FindRecurrence for an entry of the product of the Hilbert 
matrix and its inverse. 

Lemma 4.1. The summand P(n,i,m,j) satisfies the recurrence relation 

- f n -i+ifn+i-2\P(n,i- l,m,j) - P(n -l,i- l,m, j)J + 

(-l) n + i fl 1 (p(n,i,m,j)-P(n-l,i,m,j)) =0, (3) 
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and the sum p(n,i,m) satisfies the recurrence relation 

- fn-i+ifn+i-2 (p( n , i - l,m) - p(n - l,i - 1, m)\ + 

(-l) n+i fli (p(n, i, to) - p(n - 1, i, m)) = 0. (4) 

Proof. Write each of the terms in Eq. (||) as a multiple of P(n — l,i — l,m,j) to 
get the equation 

- fn-i+ifn+i-2(P(n, i - l,m,j) - P(n l,m,j)) 

{P(n, i, m, j) -P(n-l,i, to, j)) 

fn+i-2 



fn—i+lfn—j fi+j — 1 

where 



M(n,i,j)P(n-l,i-l,m,j) (5) 



M(n,i,j) = fn+i-lfn+j-l.fi+j-2 + fn-ifn-jfi+j-2 

+ ( — fn+i-ifn+j-lfi+j-l + fn-i+lfn-jfi+j-1- (6) 



It suffices to show that M(n, i,j) = 0. But this follows from the standard Fibonacci 

identities fn-ifi+j-2 + fn-i+lfi+j-1 = fn+j-1 an d fn+i-2Ji+j-l — fn+i-lfi+j-2 = 

□ 



(-l) i+i - 2 /« 



If we can establish p{n, 1,1) = 1, p(n, l,m) = if m 1, and p(n,n,n) = 1, 
then Eq. (Q) shows that p(n 7 i,rn) = 1 if i = m and p(n,i,m) = if i ^ m, for 
1 < i, to < n. 

Case 1, to): The summand P(n, 1, to, j) satisfies the recurrence 

(-l) m+1 f n -xf n +m-2P(n,l,m- l,j) - / n / n _ m+1 P(n-l,l,m-l,i) 

+ (-l) ro / n -i/n+m-iP(n, l,m, j) + f n fn- m P(n - l,i,TO, j) = 0, (7) 

and this implies a similar recurrence for p(n, 1, to). The proof of Eq. (Q) is similar 
to the proof of Lemma 4.1. The initial values of this recurrence are p{m, 1, to) and 
p(n, 1, 1). The summand P{m, l,m,j) satisfies the recurrence 

(-l) m f m f m -iP(m, l,m,j) = Gi(m, j + 1) - G x (m, j) 

where Gi(to, j) = (— l)- J_1 / J / J _iP(m, 1, to, j). Since the support of Gi is 2 < j < 
to, this equation implies that (— l) m fmfm-ip(m> 1> TO ) = 0. Therefore, when m > 1 
we get p(to, 1, to) = 0. Finally, the summand P(n, 1, 1, j) satisfies 

(-l)"/2p( n , 1, 1, j) - G 2 (n, j + 1) - G 2 (n,j), 

where Gzfaj) = (— l)- 7_1 /|P(n, 1, 1, j). In this case, the support of Gi is 1 < j < 
n, so summing over j from 1 ton gives (—l) n f%p(n, 1, 1) = —G2(n, 1) = (— l)™/ 2 , 
implying p(n, 1, 1) = 1. 

Case p(n, n, n): The summand P(n, n, n, j) satisfies the recurrence 

P{n + l,n + l,n+ l,j) - P{n, n, n,j) = G 3 (n,j + 1) - G 3 {n,j), 
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where 

Gj( „,, . ^ +2( _ ir j ((^^.(("-^ 2 

When we sum over j, the right hand side telescopes to and the left side is p(n - 
1, n + l,rt + 1) — p(n, n, n). 



This completes the proof of Theorem 2.2 



5. Reciprocal Hankel matrices based on binomial coefficients 

In this section we will prove that certain reciprocal matrices based on binomial 
coefficients have integer entries. We will give formulas for the entries of the inverses 
of these matrices. 

Let a k = ( k + 1 ). 

Theorem 5.1. Let A(n) be the n x n matrix whose (i, j)- entry is 



n + i\ I n + k\ f i + k — 1\ (i + k\ i 



Then Aij(n) is an integer, and A(n) is the inverse of the matrix R n (ak)- 

Proof. First we show that Aij (n) is an integer. We use the well known fact that if 
a is even and b is odd, then (?) is even. If i is even, then obviously Ajj(n) is an 
integer, so assume that i is odd. Now if k is also odd, then ( l ^, k ) is even, so we may 
assume that k is even. Now one of (™^) and is even. 



Theorem 5.2 below shows that that A(n) is the inverse of the matrix R n {ak)- □ 



Let bk = bk(r) be the binomial coefficient (' c+ ^ 1 ). Suppose that r a positive 
integer and r > 3. Then the inverse of R n (bk(r)) does not always have integer 
entries, but the values of n for which the inverse does have integer entries seem to 
occur periodically. Further, when the entries are not integers, the denominators are 
divisors of r. The following conjecture is true for n < 20, r < 10, and r an integer. 

Conjecture 5.1. Suppose that r is a positive integer. The inverse of the matrix 
Rn(bk(r)) has integer entries if and only if n = (mod q) or n = 1 (mod q) for 
all prime powers q that divide r. 

We do have an explicit formula for the entries of the inverse. 

Theorem 5.2. Let B(n,r) be the n x n matrix whose (i, j) -entry is 

Bij(n,r) = 

\i+k+i (n + i + r — 2 s \(n\(n + k + r — 2\ (n\i 2 Y\^Iq i+j + l 



/W\ k J\kJ rlfji + k + l' 

Then B(n,r) is the inverse of the matrix R n (bk). 

The theorem is valid if r is an indeterminate, not just if it is a positive integer. 
Also note that P>ij(n, 1) simplifies to the (i, j)-entry of the inverse of the Hilbert 
matrix, and 5y(n, 2) is equal to Aij(n). 
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Proof. Let 



i+fe+1 



n + i + r — 2\ (n 



H(n,i,m,j,k) = (-1) 

n + ft + r-2\ /n y 2 n[r ^+j + ? 1 

so that /i(n, i, m) = J2]=i Sfc=o -^( n J *i TO ; ii i s the (i, m)-entry of £?(n, r)R n (bk). 
Then H satisfies the recurrence 

n 2 (i — m + r — l)(n — z + r — l)(n + i + r — 3)H(n — l,i — l,m — l,j, k) — 

n 2 (i — to — l)(n — i + r — l)(n + i + r — 3)H(n — l,z — l,m,j, k) + 

n 2 (i — l) 2 (i — m + l)H(n — 1, i, to — 1, j, k) — 

n 2 (i — l) 2 (i — m — r + l)H(n — m,j, k) — 

(n + r - 2) 2 (i - m + r - l)(n - i + l)(n + i - l)H{n,i- l,m-l,j,k) + 

(n + r - 2) 2 {i - m - l)(n- i + l)(n + i - i)H(n,i- l,m,j,k) - 

(n + r - 2) 2 (i - l) 2 (i -m + l)H(n, i,m- l,j, k) + 

(n + r- 2) 2 (i - l) 2 (i - m - r + l)H(n,i,m,j,k) = 0. 

(8) 

The preceeding recurrence was found by FindRecurrence. The theorem will follow 
if we can establish the correct values of h(n, 1, to), h(n, n, n), and h(n, i, 1). 
Case h(n, 1,to): Maple computes h(n, 1, 1) = 1, and it computes 

#i(n,l,m,j) = 2^H(n,l,m,j,k) = (j+m+r _ 2 \ 

fe=0 r V r ) 

Now h(n,l,m) = Hi(n,l,m, j), and with Hi(n,l,m, j) as input, the function 
sumrecursion gives the recurrence (n — l)(n — 2 + m + r)h(n, 1, m) — (n + r — l)(n — 
m)h(n — 1, 1, to) = 0. Maple gives the initial value h(m, 1, m) = 0, for to > 1. 
Case h(n,n,n): Maple computes 

, £ m (-D" + ^( 2 "r' 2 )("Tr 3 )(p 

Hi(n,n,n,j) = ^ H (n, n, n, j, k) = , ra+3+r _ 2 , 

k=o r \ r ) 

Similarly to the previous case, sumrecursion gives the recurrence h(n, n, n) — (n — 
l,n — 1, n — 1) = 0, and obviously h(l, 1, 1) = 1. 

Case h(n,i, 1): We need to do something different in this case. First, we show 
that our conjectured inverse is symmetric. Let 

S(n,i,j,k) = 

* +fe +i (n + i + r -2\(n\(n + k + r-2\(n\ i 2 UZo * + 3 + 1 



(-1) 



k J \k 



r TTi=o i + k + l 



so that P>ij(n,r) = X/fc=o ^{ n ,i,j, k). Now zeil produces the recurrence S(n + 
l,i,j,k) - S(n,i,j,k) = T(n,i,j,k+ 1) - T(n,i,j,k), where 

m/ , , -(2n + r)k 2 (i + k + r -2) 

T(n,i,j,k) = - — ■ -^27 . ; 1U j— -7T S{n,i,j,k). 

(n + r - l) 2 (n — i + l)(n — k + 1) 
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This implies that Bij(n + l,r) — Bij(n,r) — T(n,i,j,j) — T(n,i,j,0). Now Maple 
tells us that T(n,i,j,j) — T(n,i,j,0) — T(n,j,i,i) + T(n, j,i,0) = 0. This means 
that Bij(n + 1, r) — Bji(n + 1, r) = Bij(n, r) — Bji(n, r). Maple tells us that 

B in (n,r) - B ni (n,r) = 

( n+i X r ~ 2 )i(n + i + r- 3)!r(2 - r)T(2 -n-i-r) (-1)* 
r(n + i- + r- 2)W(2 -n- r)Y{2 - i - r) Y{1 - i) ' 

which implies B, n (n, r) — B n i(n,r) = 0. 

Since R n (bk) and B(n, r) are symmetric, the (1, i) entry of R n (bk)B(n, r) equals 
the (i, 1) entry of B(n, r)R n (bk). The former is Y^j=i Y^k=i U{n, k), where 

U(n,i,j, k) = + ^ ^ S(n,j,i,k). 

The function zeil produces 

w • • 7 \ -(2n + r)k 2 (j + k + r -2) . 

Y(n,i,j,h) = - — — ^ — —jf — -?-—U(n,i,j,k) 

(n + r — l) z (n - j + l)(n — k + 1) 

which satisfies 

C/(n+ l,i,j,k) -U(n,i,j,k) = Y(n,i,j,k + 1) - Y(n,i,j,k). 
Thus we have 

i-l i-l 

^C/(n+ l,i,j,k) - ^U{n,i,j,k) = Y(n,i,j,i) - Y(n,i,j,0), 
fe=i fc=i 

and Maple tells us that ^™ =1 Y(n,i,j,i) — Y(n,i,j, 0) = 0. All that remains is to 

check the initial value Y^j=i Y^k=i Ji ^) = 0- Maple tells us that 

VVm- • • M = r(1 - r)r(2z - r)r(2z + r + l)r(2i + r- l) (-!)' 

2^Z^ U ^^ K > r(-r- i) 2 r(i + r + i) 2 r(i + r) (i_i)r(-»)' 

which implies that ^}=i Y^k~=i ^(*> *i Ji fc) = when i > 1. □ 

We consider reciprocal Hankcl matrices based on one more sequence of binomial 
coefficients. Let = ( fc g 3 )- 

Theorem 5.3. Let C{n) be the n x n matrix whose (i, j)- entry is 

_ >T"V i y+fc+i fn + i + 2\fn + k + l\fi + k + l\fi + k\ i(j - k) 
ClAn) -to {i + k + l){i + k+l){ i ){ i ) 3 • 

T/ien Cij(n) is an integer, and C(n) is the inverse of the matrix R n {ck)- 

Proof. First we show that each summand of the sum which defines each entry is 
an integer. It is well known that if a = (mod 3), 6=1 (mod 3), c = 2 (mod 3), 
then (^), ("), and (*)3 are all divisible by 3. Using this fact, we find that one of the 
terms ( t+k j +1 ), f^)' or * is divisible by 3 unless i = 1 (mod 3) and k = (mod 3). 
But now n + i + 2 = n (mod 3), n + k + 1 = n + 1 (mod 3), and i + k + 1 = 2 
(mod 3). Thus 3 divides one of the terms (™+ fe l +i) or (7+^). 
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The proof that C(n) is the inverse of R n (ck) is similar to the proof of Theorem 

H Let Z(n,i,mJ,k) = {-l^K^^l)^ 1 )^^^' 80 that 

z(n,i,m) = Y^j=i Sfe=o Z(n, i, m, j, k) is the (i,m) entry of C(n)R n (ck)- Then Z 
satisfies the recurrence 

(n — i + + i + 1) ( ^(ti — 1, i — 1, m, j, k) — Z(n, i — 1, m, j, fc)J + 

i(i - l)f-^(n - l,i,m,j,k) - Z(n,i,m,j, kf) 



0. 



Now the proof proceeds similarly to the proof of Theorem 5.2, except that we don't 



have to do the difficult initial value m = 1. □ 

One might wonder whether there isn't a simpler formula than the one we give 
for Bij (n, r). If we fix i and j and consider By as a polynomial of n, then it usually 
has an irreducible factor of degree min{2i — 2, 2j — 2}. Thus it seems unlikely that 
one could avoid the sum in the given formula. The next section suggests that the 
given sum is the 'right' way to describe Bij(n, r). 

6. Reciprocal Hankel matrices based on Fibonomial coefficients 

Remarkably, by changing the exponent of —1 and changing the binomial co- 
effiecients to Fibonomial coefficients in the formula for Bij, we get a formula for 
the entries of the inverses of reciprocal Hankel matrices based on Fibonomial coef- 
ficients. 

Let dk = dk{r) be the Fibonomial coefficient ^( fe+ r _1 ))- 
Conjecture 6.1. Let D(n,r) be the n x n matrix whose (i, j)- entry is 
3-1 

Dij=Dij(n,r)=^2(-ir^ 

k=0 




k ) J F r n[=0 F i+h+l 



Then the D(n,r) is the inverse of the matrix R n (dk) 



We have verified this conjecture for n < 16 and r < 10. (We assume that r is 
a positive integer.) We also observe that the inverse of a reciprocal Hankel matrix 
based on Fibonomial coefficients has integer entries exactly when the corresponding 
reciprocal Hankel matrix based on binomial coefficients has integer entries. This 
may just be a consequence of known divisibility properties of the Fibonomials. 
It seems likely that this conjecture may be proved by combining the methods of 



the proofs of Theorem 2.2 and Theorem 5.2, and that it may be extended to the 



corresponding sequence of cc-Fibonomial coefficients. 
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